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Abstract 

A  numerical  solution  is  presented  for  the  natural  convection  heat  transfer  from  two  vertical  fins  using  a  spectral 
finite  difference  method.  Virtual  distant  boundary  conditions  for  two  bodies  that  are  compatible  with  plume 
behavior  and  with  an  overall  continuity  condition  are  introduced.  A  boundary-fitted  coordinate  system  is  formed. 
Streamlines,  isotherms,  mean  Nusselt  numbers  and  drag  &  lift  coefficients  are  presented  for  a  variety  of 
dimensionless  parameters  such  as  a  Grashof  number  and  a  Prandtl  number  at  a  steady-state.  Extensive 
effectiveness  of  a  spectral  finite  difference  method  was  established. 
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1.  Introduction 

Laminar  natural  convection  has  many  engineering  applications.  In  particular,  heat  transfer  from  vertical  fins  is 
one  of  the  major  concerns  on  heat  exchangers.  There  are  a  great  amount  of  analytical  and  experimental  articles 
reported  on  natural  convection  heat  transfer  problem;  for  example,  the  flow  around  a  circular  cylinder,  a  sphere,  a 
vertical  or  horizontal  plate.  Almost  all  the  work  has  dealt  with  a  simply-connected  or  a  doubly-connected  region. 
For  most  external  natural  convection  heat  transfer  problems,  the  distant  condition  plays  an  important  role.  There 
are  many  articles  reported  on  natural  convection  heat  transfer  around  one  body  accounting  plume  behavior  [2], 
[3].  Fujii  [3]  reported  on  numerical  analysis  of  the  natural  convection  plume  above  a  horizontal  line  and  point  heat 
source  in  1963.  Mochimam  [3]  introduced  virtual  distant  boundary  conditions  that  are  compatible  with  plume 
behavior  above  uniformly  heated  horizontal  circular  cylinder  and  with  an  overall  continuity  condition. 

An  analytical  object  is  two  vertical  fins  located  in  a  row.  The  surfaces  of  two  fins  are  assumed  to  be  maintained 
at  a  uniform  high  temperature.  At  this  stage,  it  is  necessary  to  introduce  virtual  distant  boundary  conditions.  The 
virtual  distant  boundary  conditions  for  two  bodies,  the  virtual  boundary  conditions  considering  an  adjacent  body, 
and  extensive  effectiveness  of  a  spectral  finite  difference  method  [1]  are  presented. 

2.  Analysis 

Fluids  are  assumed  to  be  incompressible  and  Newtonian.  The  flow  and  the  thermal  field  are  assumed  to  be 
two-dimensional  and  laminar.  The  coordinate  system  is  nondimensionalized  with  respect  to  a  reference-body 

length  L  .  Dimensionless  variables  of  time  t  ,  a  temperature  T ,  a  stream  function  If  * ,  a  vorticity  if  are 
defined  as  follows; 


/•  evict  r-r;  w-  r_pc 

i'  V  ’  Tj-T .  vlctt  ^  vlct 

where  Gr  is  a  Grashof  number;  /,  coordinate;  V ,  kinematic  viscosity;  T  a,  an  ambient  temperature;  T  w,  a  wall 
temperature. 
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Grashof  number  Gr  is  defined  as  follows; 


Cr=L-w 

v 

where  g  is  gravitational  acceleration,  and  a  is  a  coefficient  of  thermal  expansion. 

Under  a  Boussinesq  approximation,  the  governing  equations  (a  vorticity  transport  equation,  an  energy  equation, 
the  relation  between  vorticity  and  stream  function)  can  be  written  in  non-dimensional  forms  as  follows: 


jdg  |  d(g,y/)  _  1  f  d2g  |  d2g]  |  3(7, y) 
dt  d(a,j3 )  jGr\da2  dj32  J  3(a,/?) 

J3r[3(r,^)_  i  fa2r  |  d2r] 
dt  d(a,/3)  Pr^Gri^a2  dj32  J 


,  av_n 

Jg  +  — - +  r.  —  0 

dfi- 


j  _  d(x,  y) 

d(a,j8) 


(1) 

(2) 

(3) 


where  J  is  Jacobian,  and  Pr  is  a  Prandtl  number.  Viscous  dissipation  has  been  neglected  in  the  energy  equation. 
Equations  have  been  considered  conformal  mapping  from  the  physical  domain  in  the  (x,  y )  plane  to  the 
analytical  domain  in  the  (a,  j3)  plane.  The  mapping  function  is  given  by 


x  +  iy 


K'f z  |  dn(%  +  iJ),k)cn(£  +  i7),k)  | 
n  \  sn (<^  +  iJj,k) 


ntjg  +  iri)' 
2K'K  , 


+  irj)  =  a  +  ij3 


(4) 


where  k= 0.1434,  K  and  K’  are  complete  elliptic  integrals  of  the  first  kind  with  the  complementary  parameter 


Va -k2) 


and  ^J(\  —  k'2) 


respectively.  Z  is  a  zeta  function  of  Jacobi. 


The  coordinates  and  the  geometry  are  defined  as  shown  in  Fig.  1,  where  L  is  the  fin  length.  Since  the  flow  field 
and  the  thermal  field  are  assumed  to  be  symmetric  with  respect  to  the  vertical  y-axis,  only  the  left  half  of  the 
calculate  object  is  considered  in  this  paper.  Let  dimensionless  fin  length  L  be  a  unity. 
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2.1  boundary  conditions 

Boundary  conditions  at  the  surface  of  the  left  fm  are  specified  as  follows: 

(5) 

(6) 


i/s(a0,j8)  =  ^-i/f(a0,jB)  =  0 
da 

T(a0,j3)  =  l. 


A  part  of  non-slip  boundary  condition  can  be  expressed  in  an  implicit  way  by 


~^y/(a0-d,/3)  +  JC(a0-d,/3)  =  0  (7) 

d 

where  a0  —  d  is  the  coordinate  of  the  grid  point  adjacent  to  the  surface. 

Virtual  distant  boundary  conditions  at  a  steady-state  are  introduced.  Since  the  excess  of  a  vorticity  and 
temperature  is  remarkable  in  the  plume,  the  behavior  of  a  vorticity  and  temperature,  a  stream  function  can  be 
expressed  as 


£  =  Grin,R  I  cos#0  —  cos  0 1  (/?sin  £?0  +  S )  2/5 
T„  =  Gr~u  5  (R  sin  0Q  +  S)~v  5  A4  h(A^) 
y/«,  =  Gr_3/1°  (6?  sin  <90  +^)3/5  sin  0oAf(A^)H{0) 
Go  =-Grino (RsinO0  +  S)~us  A3/^) 


(8) 

(9) 

(10) 
(11) 


cosl.4#n  +cotl.2/rsinl.4#n  .  cos0.6#n -cotl. 2;rsin0.66in  .  ,  ,  „ 

H(6)  = - 5 - -sin  0.60  + - 5 - -sin  1.40,  (12) 

sin20o  sin  26^0 

where  the  subscript  °°  means  the  virtual  distant  boundary,  S  is  a  suitable  parameter  to  be  determined.  And /  li 
are  given  by 


f”  +  -  ff"-- f'2+h  =  0 

5  5 

(13) 

h'  +  ^Prf 7i  =  0, 

(14) 

with  /?(())  =  1,  /( 0)  =  /"(())  =  =  f\°°)  =  0 .  The  prime  denotes  differentiation  with  respect  to  Ad,  . 

The  parameter  A  is  related  with  the  mean  Nusselt  number  Num  as 


A5prr  fXZ)h{Z)d£  =  Num.  (15) 

J  0 

At  the  upper  right  side  of  the  left  fin  (from  the  point  A  to  the  point  B  in  Fig.  1 .).  H  is  determined  so  as  to  be  a 
unity.  And,  the  lower  right  side  of  the  left  one  (from  the  point  C  to  the  point  D)  is  expressed  through  supposition 
of  a  uniform  flow.  The  stream  function  near  to  the  y-axis  (from  the  point  B  to  the  point  C)  is  constant  considering 
existence  of  the  right  side  fin.  On  this  virtual  boundary,  heat  flux  is  specified  to  be  zero,  and  the  vorticity  is 
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vanished.  Zero  heat  flux  boundary  conditions  are  considered  at  the  virtual  boundary  from  C  to  D  too. 


H(ff)  =  1 

at 

AB 

(16) 

y/  =  Gr-vl0(R  sin  0O 

+  S)315  sin  60Af  (A£) 

at 

~BC 

(17) 

at 

BC 

and  CD 

(18) 

O 

II 

Nn 

at 

BC 

and  CD 

(19) 

ys  =  Gr~mo(R  sin  0Q 

r..  3/5  cos#-cos#n 

+ sy  - - 

sin  0oAf(A%) 

at 

CD 

(20) 

cos#,  —  cos  Oq 


where  E,  is  a  value  of  independent  variable  when  an  argument  0  =  #, . 

At  a  partial  Neumann-type  virtual  thermal  boundary  conditions  Eqs.  (18),  a  following  equation  can  be  used. 

T(aj,j3,t)=  1  b1(aj,jB,t)  +  cT(aj  -d,j3,t)-(l  +  c)T(aj  -2 d,j3,t)}  (21) 

da  (2  +  c)d 

where  d  is  a  grid  size  near  virtual  boundary,  and  c  is  a  constant  chosen  arbitrarily[4]. 

2.2  Numerical  procedure 

Let  l// ,  C  ,T  be  expressed  as 


¥  =  Yj  Yck  cos  kp  +  Yj  ¥sk  sin  k/J 

k= 0  k= 1 

(22) 

C  =  X  £ck  cos  kfi  +  Y  Csk  sin  kP 

k=0  k= 1 

(23) 

r  =  YTck  cos kp  +  X Tsk  sin kf3  . 

k= 0  k= 1 

(24) 

Substituting  Eqs.  (22)  -  (24)  for  governing  equations,  and  decomposing  the  equations  into  cosine  terms  and  sine 
terms  gives  the  forms  as 

X  fck  («>  0  cos  kp  +  Yj  fsk  ( a ,  t)  sin  kj3  =  0 

k= 0  k=0 


which  reduces  to 


fck(a,t)  =  0,  fsk(a,t)=  0,  k>  0. 

Then  each  Fourier  component  in  f3  can  be  separated  from  Eqs.  (1)  -  (3).  Hereafter  subscripts  ck,  sk  mean  a 
wave  number  of  cosine  term,  and  a  wave  number  of  sine  term  respectively. 

All  second-order  derivatives  with  respect  to  a  are  replaced  by  three-point  central  difference  approximation. 
First-order  derivatives  with  respect  to  a  near  boundary  surface  are  used  three-point  forward  difference 
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formulae.  Time  integrations  are  replaced  by  semi-implicit  method.  With  these,  the  differential  equations  are 
converted  into  a  set  of  algebraic  equations.  These  equations  are  then  solved  iteratively  per  wave  number. 

Time  difference  approximation; 


dt  1  At 


(25) 


Spatial  difference  approximation; 


fk-idk 


d/ 

a=ak  dk_i(dk  +dk_1) 

37  2/w 


-  + 


dcr 


dk-i(dk  +dk_  j) 


fk(dk  dk_\)  fk+\dk-\ 
dk_xdk  dk(dk+dk_l) 

2fk  [  fk+ 1 

dk-\dk  dk(dk  +dk_  j) 


(26) 

(27) 


To  support  high  Grashof  number  cases,  non-uniform  grid  spacing  is  used. 


a, 


a{)  -  c 


[ sinh(k  -1  )y 
sinh  y 


+  1 


(28) 


where  OCk  is  a  coordinate  at  the  A-th  grid  point  counted  from  the  wall  along  CC  coordinate,  y  is  a  suitable 
constant  and  c  is  a  constant  to  be  determined. 


2.3  Dimensionless  force 

Dimensionless  force  acting  on  a  fin  surface  along  positive  x  direction  is  defined  as  a  drag  coefficient  Ca.  And, 
dimensionless  force  acting  on  a  fin  surface  along  positive  y  direction  is  defined  as  a  lift  coefficient  CL 


Cr 


C, 


dx  dC 

- +  —  y 

djB  da  - 


dp 


-!  f»|  oj_ 

•Ur  ^  7*R 


dy  d£ 


dp 


4Gr]-\  dp  da'  j 

A  mean  Nusselt  number  Num  at  the  surface  of  the  left  heat  source  is  defined  as  follows: 
Nu,„  =  -n- 


da 


(29) 

(30) 


(31) 


2.4  Convergence  criteria 

To  seek  a  steady-state  solution,  the  iterative  procedure  is  terminated  when  the  following  convergence  criteria  Eqs. 
(32),  (33)  are  satisfied  at  the  surface  of  the  fin. 

<  10-6  (32) 

<  10-6  (33) 


C  -C 

^D(new)  ^D(old) 

r 

^D(new) 

C  -C 

'“'L(  new)  ''''L(old) 

c 

^L(new) 
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Nu 


m(new) 


■Nu 


m(old) 


Nu 


m(new) 


<10 


-6 


(34) 


where  the  subscripts  (new)  and  (old)  denote  the  current  value  and  the  value  in  the  past  respectively. 

The  number  of  partitions  along  CX  ,  and  tine  wave  number  to  seek  a  steady-state  solution  is  determined  so  that 
the  variation  of  dimensionless  forces  does  not  exceed  one  percent.  In  this  paper,  the  number  of  partitions  along  CX  , 
and  the  wave  number  are  used  in  the  range  of  100  to  200,  and  in  the  range  of  30  to  40,  respectively. 

3.  Results  and  discussion 

A  numerical  solution  for  steady-state  natural  convection  heat  transfer  accounting  plume  behavior  around  two 
vertical  fins  whose  temperatures  are  maintained  at  a  higher  temperature  than  that  of  the  ambient  has  been 
obtained,  against  a  Rayleigh  number  in  the  range  of  0.7  to  70000  for  air  (Pr=  0.7)  and  water  (Pr= 7). 


(a)  (b) 

Fig.  2.  (a)  Steady-state  streamlines;  Ay/  =0.24,  (b)  isotherms;  AJ=0.1  at  Pr  =0.7,  Gr  =1 


Fig.  3.  (a)  Steady-state  streamlines;  Ay/  =0.07,  (b)  isotherms;  AT  =0.1  at  Pr  =0.7,  Gr  =1000 
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Fig.  4.  (a)  Steady-state  streamlines;  Ay/  =  0.135,  (b)  isotherms;  AT  =0.1  at  Pr  =7,  Gr  =1 


Fig.  5.  (a)  Steady-state  streamlines;  Ay/  =  0.035,  (b)  isotherms;  AT  =0.1  at  Pr  =7,  Gr  =1000 


Fig.  6.  (a)  A  lift  coefficient  CL,  (b)  a  drag  coefficient  CD  against  PrGr 
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Fig.  7.  A  mean  Nusselt  number  Nu„  against  PrGr 


Steady-state  streamlines  and  isotherms  are  shown  in  Fig.  2-5.  Thermal  fields  similar  to  solid  heat  conduction 
are  obtained  in  case  of  a  low  Grashof  number  (see  Fig.  2  (b)).  A  steady-state  buoyancy  plume  above  the  fin  is 
formed  as  an  increase  of  a  Grashof  number. 

The  lift  coefficient  and  the  drag  coefficient  at  a  steady-state  for  various  Grashof  numbers  and  Prandtl  numbers 
are  presented  in  Fig.  6.  The  mean  Nusselt  number  of  the  present  at  a  steady-state  is  proportional  to  a  Grashof 
number  to  the  power  0.2(see  Fig.  7).  Clearly,  the  mean  Nusselt  number  of  water  always  indicates  higher  heat 
transfer  rate  than  that  of  air,  and  the  lift  coefficient  of  water  within  the  results  of  this  investigation  is  bigger  than 
that  of  air  in  case  of  a  same  Rayleigh  number.  However,  the  drag  coefficient  of  water  is  smaller  than  that  of  air  in 
case  of  a  same  Rayleigh  number.  The  mean  Nusselt  number  and  the  lift  coefficient  are  in  good  agreement  with 
the  solution  of  the  boundary  layer  equations. 


4.  Conclusion 


Steady-state  laminar  natural  convection  heat  transfer  around  two  vertical  fins  using  the  virtual  distant  boundary 
conditions  and  the  virtual  boundary  conditions  was  studied  numerically  by  a  spectral  finite  difference  method  for 
various  Grashof  numbers  and  Prandtl  numbers.  The  virtual  distant  boundary  conditions  and  the  virtual  boundary 
conditions  for  two  bodies  could  demonstrate  plume  behavior  very  well.  Extensive  effectiveness  of  a  spectral  finite 
difference  method  was  established. 
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